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The band structure of two-dimensional artificial superlattices in the presence of (Rashba-type) 
spin-orbit interaction (SOI) is presented. The position and shape of the energy bands in these 
spintronic crystals depend on the geometry as well as the strength of the SOI, which can be tuned 
by external gate voltages. For finite mesoscopic arrays we show that their conductance properties 
can be understood from these spin-dependent band diagrams. 
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I. INTRODUCTION 


Infinite periodic structures - although in the strict 
sense they do not exist - provide several fundamental 
concepts that determine the physical properties of var¬ 
ious systems. E.g., the notion of bands is essential in 
the theory of solids. The concept of ’’crystals” has also 
been used e.g. in the context of photonic bandgap mate¬ 
rials 0 , 0 - In that case artificial periodic structures have 
been fabricated with a photonic band structure. Here 
we show how both the energy and spin of electrons can 
be engineered leading to the concept of spintronic crys¬ 
tals. The band structure of these crystals is flexible in 
the sense that instead of being completely determined by 
the geometry, it is tunable by gate voltages. 


We consider two-dimensional superlattices, which can 
be fabricated from e.g. InAlAs/InGaAs based het¬ 
erostructures (3) or HgTe/HgCdTe quantum wells itj, 
where Rashba-type (Hj spin-orbit interaction (SOI) is 


present. This effect, which is essentially the same as 
the one which causes the fine structure in atomic spec¬ 
tra, results in spin-dependent transport phenomena. Ex¬ 
periments have demonstrated that the strength of this 
type of SOI can be controlled by external gate voltages 
[q, [?] . Small periodic structures like 5x5 ring arrays 
have recently been realized experimentally {8| and have 
been described theoretically 19 11101 as well. Finite chains 
of quantumrings 11 ], ladder [ 121 ] and diamond-like ele¬ 
ments [l3|, E3] have also been investigated, as well as ar¬ 
tificial crystal-like structures Em in contexts different 
from the current one. The spin transformation proper¬ 
ties of finite networks su gge st various possible spintronic 
17] applications as well 13, 18. E]. 


Conductance calculations show - already in relatively 
small finite lattices - signatures of a possible underlying 
band structure, as one can identify SOI-dependent en¬ 
ergy regions where the system is completely opaque for 
the electrons. In the present paper we show that these 
” non-conducting stripes” m are directly related to the 
relevant band gaps. The general two-dimensional (2D) 
lattice that we will consider is shown in Fig. [1} We found 
that the details of the band structure strongly depends 


on the geometry of the unit cell, but the overall scaling 
properties, as well as the remarkable SOI dependence are 
general. 

The characteristic energies in our problem are much 
smaller than e.g. the usual electronic band gaps in semi¬ 
conductors. This is essentially due to the differences in 
the lattice constants, typical ’’bond lengths” have to be 
compared to spatial periodicity of the artificial superlat¬ 
tices, which is around 10 nrn. Thus the nanometer-scale 
translational symmetry induced subbands that we de¬ 
scribe in the current paper can be considered as a fine 
structure of the usual electronic bands, which becomes 
important at low temperatures. 




FIG. 1: Two-dimensional parallelogram lattice (a) and a finite 
array (b), where the relevant parameters (that may lead to 
square or rectangular lattices as well) and some possible input 
and output leads are shown. Electrons move along the lines 
connecting the junctions (full circles). 


II. CALCULATION OF THE SPIN-DEPENDENT 
BAND STRUCTURES 


The Hamiltonian of a narrow quantum wire in the x- 
y p lane with Rashba-type SOI can generally be written 
201 as 
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where the unit vector n points to the chosen positive 
direction along the wire, and we introduced the charac¬ 
teristic kinetic energy Ml = H 2 /2m*a 2 (with a being one 
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of the lattice constants, see Fig. [TJ. The strength of the 
SOI is given by ui = a/a, where the Rashba parameter 
a i s tunable by gate voltages, and s denotes the (di¬ 
mensionless) length variable along the wire measured in 
units of a. 

Independently from the direction of the wire, the en¬ 
ergy levels of H form a continuum, and the spin direction 
of the eigenspinors depend on the SOI strength. A given 
energy eigenvalue is fourfold degenerate due to the two 
possible propagation and spin directions. The structure 
shown in Fig. [T) corresponds to two Hamiltonians with 
an angle of 7 between the corresponding n vectors. In 
order to find an eigenstate for the whole geometry, the 
solutions have to be fitted at the junctions. We imply 
Griffith’s boundary conditions 21], that is, the net spin 
density current at the junction has to be zero and we also 
require the continuity of the spinor valued eigenfunctions. 

In the case of an infinite periodic structure, we look 
for Bloch-wave solutions [exp(ik ■ r)^p(r) with lattice- 
periodic spinors </?)], which means an additional, special 
boundary condition. The consequence is an energy spec¬ 
trum with a specific structure: e.g., there will be no so¬ 
lutions in certain energy ranges. Determination of the 
band structure means finding triples {.E(k), k±, that 
correspond to a Bloch-wave eigenspinor of the problem. 
As usually, we find that -E(k) is a multi-valued function 
of the two dimensional wave vector k, and we can identify 
infinitely many surfaces in this function. These surfaces 
(bands) do not overlap unless in the presence of a sym¬ 
metry induced degeneracy. 

Note that the model above assumes single mode prop¬ 
agation, which is a reasonable approximation for narrow 
conducting wires. Taking the finite width of these chan¬ 
nels into account leads to qualitatively the same results, 
with considerably increased computational costs. Addi¬ 
tionally, subbands related to the transversal modes have 
already been analyzed in detail (see e.g. Ref. |22j), thus 
using the current model we can focus on the band struc¬ 
ture induced by the periodicity of the lattice. 


III. BAND SCHEME AND CONDUCTANCE 
PROPERTIES 

For a given two-dimensional wave vector k the energy 
eigenvalues can be written as E n ^ m { k), where the band 
indices n and m are related to the spatial periodicity of 
the exp(k • r) waves in the unit cell along the two lat¬ 
tice directions. These energies scale essentially with the 
square of n/a and m/b. The same phase relations at the 
boundaries can hold with e.g. n and n+l waves along the 
direction of one of the lattice vectors in the unit cell, and 
the dominant contribution (omitting SOI corrections) of 
these solutions to the energy is proportional to n 2 /a 2 and 
(n + l) 2 /a 2 . That is, the bands have in general a double 
quasiperiodic structure. However, when a ~ b, a repeti¬ 
tion of a small number of E n rn (k) surfaces provides the 


complete band structure. 
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FIG. 2: Band structure of a lattice with 7 = n/2 and b/a = 
1.03, with strength of the SOI (a) zero and (b) ui/Sl = 5.0. 
Cross sections at = 1.0 are shown for several additional 
values of the SOI strengths in panel c). The thin dotted 
lines guide the eye by showing the energy range for which 
the bands are ’’essentially the same,” i.e., they continuously 
transform into each other when the SOI strength is changed. 
Note that the levels between the dotted lines correspond to 
the two dimensional plots in panels a) and b). 


Figure [3(a) shows four bands for a rectangular lat¬ 
tice with b/a = 1.03 and zero SOI. (Note that energy 
is measured in units of HQ, which, for a = lOnm in In- 
AlAs/InGaAs based heterostructures is of the order of 
meV). The four bands seen in this figure are quasiperi- 
odically repeated. The lattices shown in Fig. [1] have unit 
cells with four distinct boundary points, that is, four 
leads connect them to the neighboring cells. The oppo¬ 
sitely situated boundary points are equivalent in a crys¬ 
tal, thus any measurable physical quantity has to have 
the same value at these points. Particularly, the currents 
carried by the opposite leads should be the same. That 
is, the sign of the currents at the four leads can be written 

schematically as + + ++, H- 1 —, —I-b and- 

(where the leads that correspond to the ± signs follow 
each other in a clockwise order.) The four bands seen in 
Fig. (Ha) correspond to these four possible current con¬ 
figurations. For nonzero SOI, all these bands split into 
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two due to the spin dependence of the interaction and 
as it is shown in Fig. [2J the strength of the SOI modi¬ 
fies considerably both the position and the width of the 
allowed/forbidden bands. 

In Fig. 0(c) cross sections of the band structure are 
plotted for different values of the SOI which clearly shows 
the gradual splitting of the levels as the SOI gets stronger. 
Additionally, when we identify the bands that continu¬ 
ously evolve from/into each other when the strength of 
the SOI [characterized by the parameter w in Eq. (QJ] is 
changed, we notice an overall decrease of the energies [see 
the dotted lines in Fig. [2][c)]. This is due to the SOI in¬ 
duced splitting of the lowest band, resulting in a decrease 
of the lowest possible energy when u increases. 



FIG. 3: The effect of the lattice geometry on the band struc¬ 
ture for lj/Q — 1 . 0 , k\a = 1 . 0 . 


Figure |U[a) shows a contour plot of the conductance as 
a function of the energy and the SOI strength for a rect¬ 
angular 15 x 15 array. We clearly notice the appearance 
of stripes (the position and width of which depend on the 
SOI strength) of zero conductance. In these regions the 
array is completely opaque for the electrons. Addition¬ 
ally, for the 15 x 15 array, these ’’non-conducting stripes” 
[HI coincide with the band gaps obtained from a calcu¬ 
lation assuming an infinite structure with the same local 
geometry. In order to visualize this fact, we projected 
the band structure on the energy axis to obtain the lim¬ 
its between allowed and forbidden energy regions [see the 
grey areas in Fig. [4] (b)[. Already for a 3 x 3 network, we 
can see some signatures in G(E) of the band structure, 
but for a 7 x 7 array the positions of the zero conductance 
energy ranges are practically the same as the band gaps. 
Having introduced random scattering centers as it is dis¬ 
cussed in Refs. mm , we observed that the widths of 
the band gaps decrease only 10%, even when dephasing 
is so strong, that 100% degree of spinpolarization (input) 
drops to 20% (outputs). 



The band structure strongly depends on the underlying 
geometry as illustrated in Fig. [3] This emphasizes that 
besides the tunability of the band structure by external 
gate voltages (that modify the strength of the SOI), ge¬ 
ometry is also an important additional degree of freedom. 

Results based on infinite structures usually have im¬ 
plications also on large, but finite systems. Now we cal¬ 
culate the conductance of arrays consisting of N x N 
unit cells [see Fig.QJb)] using the Landauer-Biittiker 22] 
formula: 




where the sum runs over the possible outputs. T?{E) 
[T™(E)] refers to the transmission probability at the rel¬ 
evant output for spin-up (spin-down) input in the chosen 
quantization direction. These probabilities are calculated 
by solving the eigenvalue problem for the whole network 
at a given energy E imposing the appropriate boundary 
conditions, e.g. at the input we have a spin-up (or spin 
down) incoming wave and a possible reflected one, while 
at the outputs only outgoing waves appear. 



FIG. 4: Top: Contour plot of the conductance measured in 
units of Go = e 2 / h for a 15 x 15 array with b/a = 2,7 = 7 r /2 as 
a function of the energy and the SOI strength. Panel b) shows 
the conductance of 3 x 3, 7 x 7, 15 x 15 arrays where the light 
grey shading indicates the energy gaps in the corresponding 
infinite superlattice. 
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The most straightforward possible application is en¬ 
ergy filtering: E.g., when a s=s b there are very narrow 
bands the width and position of which can be controlled 
by the SOI strength. (For InAlAs/InGaAs based het¬ 
erostructures and a 9 x 9 network with a = 10 nm, 
b/a = 1.005, an energy range around 5 meV is trans¬ 
mitted in the middle of a 100 meV wide band gap.) 

Additionally, as the width of all the band gaps can be 
controlled simultaneously by the SOI strength [see the 
almond-shaped minima in Fig. [4] (a)], an important high 
temperature effect can be foreseen: even when the input 
has a broad energy distribution, conductance is modu¬ 
lated by the SOI. For the same device mentioned above, 
the conductance changes 20 % of its average value when 
the SOI strength is varied in experimentally achievable 
range. For a 13 x 13 network the modulation is around 
40%. (Note that in the framework of our model, at ’’high 
temperatures” transversal modes other than the ground 
state in this direction should not be excited. However, 
the physical reasons of the result above are valid also for 
multimode propagation.) 

For non-square lattices the geometrical anisotropy 
leads to anisotropy in the band structure (see Fig. [3]), 
and consequently also in the conductance properties. For 
a 9 x 9 lattice with b/a = 2, 7 = 7 r/ 4 , the difference of 
the transmission probabilities in the x and y directions - 
depending on the SOI strength can be zero, or as large 
as ±0.8, so that the higher one is above 0.95. 

Besides the SOI controlled phenomena discussed 
above, finite arrays can also perform various spin trans¬ 
formations. Apart from spin rotations that can also be 
done with smaller devices, the arrays considered here are 
also versatile spintronic devices: e.g., according to our 
calculations, the network described in the previous para¬ 
graph can deliver oppositely spinpolarized outputs from 
a completely unpolarized input, when the output leads 
are situated at the middle of the sides of the network [see 
Fig. mb)]. Note that these properties are similar to that 
of a ring arrays [HI, and although in the current case the 
transmission probabilities are lower than unity (but still 
around 50%), the point that there is no need for the local 
modulation of the SOI strength makes the arrays consid¬ 
ered in this paper more promising from the viewpoint of 
the possible applications. 

IV. SUMMARY AND CONCLUSIONS 

In this paper we investigated two-dimensional super¬ 
lattices in which the propagation of the electrons is deter¬ 
mined by the interplay of the geometry and the spin-orbit 
interaction (SOI). We calculated the band structure of 
these artificial crystals, and showed that by changing the 
SOI strength in the experimentally achievable range, the 
band scheme can be modified qualitatively, e.g., forbid¬ 
den energy ranges can become allowed and vice versa. 
Comparing the band structure with the conductance 


properties of finite systems, we found that already for 
relatively small arrays, forbidden bands are clearly seen 
in the conductance. This effect - being based on robust, 
bulk-like properties - can be useful for the development 
of spin-dependent electronic devices. Several possible ap¬ 
plications were given, including the strong modulation of 
the conductance at moderate temperatures. 

Acknowledgments We thank M. G. Benedict and 
F. Bartha for useful discussions. This work was sup¬ 
ported by the Flemish Science Foundation (FWO-V1), 
the Belgian Science Policy (IAP) and the Hungar¬ 
ian Scientific Research Fund (OTKA) under Contracts 
Nos. T81364, M045596. P.F. was supported by a 
J. Bolyai grant of the Hungarian Academy of Sciences. 


* Electronic address: foldi@physx.u-szeged.hu 

[1] E. Yablonovitch, Phys. Rev. Lett. 58, 2059 (1987). 

[2] S. John, Phys. Rev. Lett. 58, 2486 (1987). 

[3] T. Koga, J. Nitta, T. Akazaki, and H. Takayanagi, Phys. 
Rev. Lett. 89 , 046801 (2002). 

[4] M. Konig, A. Tschetschetkin, E. M. Hankiewicz, 
J. Sinova, V. Hock, V. Daumer, M. Schafer, C. R. Becker, 
H. Buhmann, and L. W. Molenkamp, Phys. Rev. Lett. 
96 , 076804 (2006). 

[5] E. I. Rashba, Sov. Phys. Solid State 2 , 1109 (1960). 

[6] J. Nitta, T. Akazaki, H. Takayanagi, and T. Enoki, Phys. 
Rev. Lett. 78 , 1335 (1997). 

[7] D. Grundler, Phys. Rev. Lett. 84, 6074 (2000). 

[8] T. Bergsten, T. Kobayashi, Y. Sekine, and J. Nitta, Phys. 
Rev. Lett. 97 , 196803 (2006). 

[9] Z. Zhu, Y. Wang, K. Xia, X. C. Xie, and Z. Ma, Phys. 
Rev. B 76, 125311 (2007). 

[10] O. Kalman, P. Foldi, M. G. Benedict, and F. M. Peeters, 
Phys. Rev. B 78 , 125306 (2008). 

[11] B. Molnar, P. Vasilopoulos, and F. M. Peeters, Phys. 
Rev. B 72 , 075330 (pages 7) (2005). 

[121 H.-X. Wang, S.-J. Xiong, and S. Evangelou, Phys. Lett. 
A 356, 376381 (2006). 

[13] A. Aharony, O. Entin-Wohlman, Y. Tokura, and S. Kat- 
sumoto, Phys. Rev. B 78, 125328 (2008). 

[14] S. Sil, S. K. Maiti, and A. Chakrabarti, Phys. Rev. B 79 , 
193309 (2009). 

[15] M. Onoda and T. Ochiai, Phys. Rev. Lett. 103, 033903 
(2009). 

[16] S. Bcllucci and P. Onorato, Eur. Phys. J. B 73, 215 

( 2010 ). 

[17] D. D. Awschalom, D. Loss, and N. Samarth, Semicon¬ 
ductor Spintronics and Quantum Computation (Springer, 
Berlin, 2002). 

[18] P. Foldi, O. Kalman, M. G. Benedict, and F. M. Peeters, 
Nano Lett. 8, 2556 (2008). 

[19] O. Kalman, T. Kiss, and P. Foldi, Phys. Rev. B 80, 
035327 (2009). 

[20] F. E. Meijer, A. F. Morpurgo, and T. M. Klapwijk, Phys. 
Rev. B 66, 033107 (2002). 

[21] S. Griffith, Trans. Faraday Soc. 49 , 345 (1953). 

[22] S. Datta, Electronic transport in mesoscopic systems 
(Cambridge University Press, Cambridge, 1995). 



5 


[23] P. Foldi, O. Kalman, and F. M. Peeters, Phys. Rev. B 
80, 125324 (2009). 



